On the gravitational instability of a dissipative medium. 
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This paper shows that the ordinary Jeans wave number can be obtained as a limiting case 
of a more general approach that includes dissipative effects. Corrections to the Jeans critical mass 
associated to viscosity are established. Some possible implications of the results are finally discussed. 



I. INTRODUCTION 



Structure formation in fluid systems is closely related to the existence of growing modes while considering density 
fluctuations in the hydrodynamical regime. Relevant growing modes at very small wave numbers have long ago been 
identified in gravitational systems, leading to key concepts in present astrophysics such as the Jeans wave number and 
the Jeans mass [Q . These parameters correspond to the conditions that a density fluctuation must satisfy in order to 
grow up to the point of forming a structure. 

Entropy production plays an important role in realistic descriptions of thermodynamical processes, and some 
discussion of its physical sources may yield interesting results in the analysis of the time evolution of statistical 
fluctuations in astrophysical systems. This work is mainly devoted to the study of the dispersion relation that arises 
from the fluctuation analysis of the thermodynamical variables of a simple dissipative fluid under the influence of a 
gravitational field. The analysis is reviewed in the case of non- vanishing viscosity and different features of the fluid are 
established. In particular, this work shows the possibility of viscosity-affected growing modes for some astrophysical 
fluids, which may be of interest for both astrophysicists and cosmologists. 

This paper is divided as follows: section two is dedicated to the hydrodynamical model analysis and the establish- 
ment of the dispersion relation in which the main results are based. Section three is devoted to the mathematical 
treatment needed for its solution. The Jeans instability is recovered in the presence of viscosity for typical situations 
present in protogalaxics. Finally, a discussion of the results obtained and an outline of future work are included in 
section four of the paper. 



In order to establish the basic equations that lie at the core of this work, we remind the reader that the so-called 
Navier-Stokes-Fourier equations of hydrodynamics for a simple fluid arise from the structure of what is now called 
Linear Irreversible Thermodynamics (LIT) by supplementing the two conservation equations for mass and momentum, 
respectively and the balance equation for the internal energy, with additional information. Indeed, these equations 



II. HYDRODYNAMIC MODEL 
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where ^ = £ + u i J- [ . 

Here, p(x J ,t), u l (x J ,t) and e(x 3 ,t) are the local density, velocity and internal energy, respectively, S 1 - 7 the momen- 
tum current (or stress tensor) and Q 3 the heat flux. All indices run from 1 to 3. In general, for isotropic fluids, 
S JJ = p8 lJ + t 1 3 \ where p is the local hydrostatic pressure and t 13 the viscous tensor (5 lJ is the unitary dyadic). 
Notice that Eqs. (Q-||) contain fifteen unknowns (including the pressure) and there are only five equations, so the 
system is not well determined. If we arbitrarily choose to describe the states of the fluid through the set of variables 
p(x 3 ,t), u l (x 3 ,t), T(x J ,t), where T is the local temperature, we need nine dynamic equations of state (or constitutive 
equations) relating E* 3 and Q 3 to the state variables plus two local equations of state p = p(p, T) and e = e(p,T). If 
according to the tenets of LIT we choose the so-called linear constitutive laws, namely 

r ij = -2r)a ij -C05 ij (4) 



with 77 and £, the shear and bulk viscosities, respectively, a 13 the symmetrical traceless part of the velocity gradient, 
9 = and k being the thermal conductivity. 

Eqs. (|j) and (||) are the well known constitutive equations of Navier-Stokes and Fourier, respectively. Substitution 
of these equations into Eqs. (|l|-||) yields a set of second order in space, first order in time non-linear coupled set of 
partial differential equations for the chosen variables p, u l and T. This set, which the reader may seek in the literature 
P|, ||, is the so called Navier-Stokes- Fourier system of hydrodynamic equations. The non-linearities appearing in such 
equations have two sources, the inertial terms u l -g~ arising from the hydrodynamic time derivatives, plus quadratic 
terms in the gradients of velocity arising from Eqs. (Q) and (^|). Moreover, it should be mentioned that this set 
of equations is consistent with the second law of thermodynamics, the Clausius uncompensated heat , or entropy 
production, is strictly positive definite. 

Nevertheless, for the purpose of this paper, this set of equations is too complicated. In order to deal with fluctuations 
around the equilibrium state, one assures that for any of two state variables, call them X(x 3 ,t) one can write that, 

X{x\t) =X a + 5X(x\t) (6) 

where X a is the equilibrium value of X and 5X the corresponding fluctuation. Neglecting all terms of order (5X(x : ' , t)) 2 
and higher in the NSF non-linear set one finds the linearized NSF equations of hydrodynamics ||- j7j, 

^(Sp) + pS = (7) 
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since Uk — 0, = Suk and = (ffr-)p , kt = Po{^-)t and Dt = k c is the thermal diffusivity. C p and C v are 
the specific heats at constant pressure and constant volume, respectively. 

Note now that from Eq. (||), Trans = Rot(uk) uncouples from the hydrodynamic modes, whence, one formally 
arrives at the desired set, namely: 

^-(6p) + Po 9 = (10) 
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assuming that Fk — — Jpr ■ Also, D v = (jm + <;)-£-■ Now, we notice that (Sp) = —p 9 so we can reduce this set to 
only two equations, Eq. (|ll]) and Eq. ( We now perform a few minor transformations introducing the speed of 
sound of the fluid C 2 = (f^J through the relationship 

C 

k T = -fk s = —?-k s (13) 
where k s is the adiabatic compressibility ^§^J whence, fcy = This leaves us finally with the set: 

dt 2 7 dx % dxi 7 dx l dxi dx l dxi dt dx l dx^ 

where we have used that C D — CL = ^ 2 T ° . 

Eqs.( pT[]l5| ) form a set of coupled equations for the density and temperature fluctuations in the fluid under the action 
of a conservative force whose nature need not to be specified for the time being. They are the basis for studying the 
properties of the time correlation functions of thermodynamic fluctuations. Those of the density will be of particular 
interest here. 



III. SOLUTION TO THE HYDRODYNAMIC EQUATIONS 

The results derived in the previous section are far from being new. Aside from the term V 2 ^!> which arises from 
the presence of an external conservative force, they are identical to the ones that have been widely discussed in the 
literature. The question here is if the gravitational potential introduces any substantial modification in the correlation 
functions for the thermodynamic fluctuations. To examine this possibility we recall that, if we are considering only 
the fluctuations, the gravitational potential satisfies the Poisson equation: 

6ij WBxl = -^ ^ 

Now, equation ( |l5| ) reads 
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The solution to Eqs.(15) and ( |l7| ) proceeds in the standard fashion. We reduce them to a set of algebraic equations 
by taking their Laplace-Fourier transform, choose to set the static temperature fluctuations equal to zero, eliminate 

the temperature leading to an equation for 8p (k, sj , which is the ratio of two polynomials in s . To compute Sp (k, tj 

one must take the inverse Laplace transform of the former quantity, which demands the knowledge of the roots of the 
denominator, which is a cubic equation in s 3 (dispersion equation). Next, we proceed with intuition, we assume that 
the Rayleigh peak is nos substantially modified by the external potential, so that according to the standard theory, 
one of the roots of the cubic equation is s = —Dt k 2 , where Dt is the thermal diffusivity. This allows us to simplify 
the denominator to the form, 

{s + D T k 2 ){s 2 + D v k 2 s + C^k 2 - AirGpo) = (18) 



the two roots of the quadratic equation are: 

D v k 2 . 
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Iff AirGpo = and viscosity dominates over the term in k 2 , this result reduces to the one giving rise to the standard 
Brillouin peaks which correspond to density fluctuations of the type H- M, 
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which are the acoustic modes damped by the Stokes-Kirchhoff factor D v . 



(20) 



On the other hand, if 4irGp ^ 0, the threshold value for k distinguishing between damped oscillations and growing 
modes is given by 
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Eq. (j22j) is a generalization of Jeans wave number when dissipative effects due to viscosity are non-negligible, and is 
the main result of the paper. We can note that, if 47T Qi° v <§C 1, and taking the — sign for the square root, we have: 
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This is the square value of the Jeans wave number § §. It was derived by Jeans in 1902 and rederived in many 
other waves by several authors H . Here we simply show that it is almost a trivial consequence of LIT. 



Let us expand the square root in Eq. ( j22|) up to second order in order to have a better grasp of the effects of 
viscosity on the Jeans wave number. Now we have that, 
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or, using Eq. (E4): 
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so that the dimensionless number 



irG 



is a good measure of the second order corrections to the Jeans wave number due to viscous effects. 
For the corrected Jeans mass we have Hi: 
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M c = Mj(l+e)-3 
where the standard Jeans mass Mj is given by ||: 
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Mj = - 
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One should compare these results with Weinberg || and Peebles p. 108 ||. The question here is if in the primeval 
Universe, the plasma had a substantial viscosity. 

Eq. ( p0| ) deserves more attention. As shown by Weinberg Q a similar result follows from a more intuitive calculation 
addressing the behavior of the dispersion relation in sound waves. Using numerical values for the density of barions 
(protons) and the velocity of sound in a background of photons and protons, he is able to estimate the temperature 
range for which the protogalactic mass can grow to 10 11 solar masses. Similar reasoning may be applied to Eq. j29| ) 
and, of course, predict a new Mj. The question here is if the viscous corrections have any real meaning at all or they 
just turn out to be an academic curiosity. For instance, in accretion disks (^] p ~ 4 x 10~ 4 ^j-, r\ can be estimated as 

high as 10 10 Pa • s and C a « 5 x 10 3 ^ . This would yield a value of D v « 10 14 ^, so that 

2D 2 

e = -^irGpo w 10- 1 (31) 
the estimated correction is at the order of 10%. 



IV. FINAL REMARKS 



Dissipation is always present in real processes. Whenever the thermodynamical functions of a self-gravitating viscous 



uilibrium condition for the stability 
that measures how important this 



system fluctuate around an equilibrium value, dissipation will affect the general eq\ 
of those fluctuations. We have been able to establish a simple expression, Eq. 
correction will be, assuming that the Rayleigh peak remains essentially unaffected by the the gravitational field. If 
this were not the case, then a full treatment of the dispersion relation would be needed, involving interesting effects 
of the gravitational field in the Rayleigh-Brillouin spectrum, for the case of damped systems. This modifications 
may have some importance in astrophysics, since natural astrophysical masers are already well-known [|l0| and its 
scattering in fairly dense systems will provide information about the transport properties of the scatterer. 
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